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Abstract 

We study the production of cc pairs and dimuons from hard collisions in 
nuclear reactions within the covariant transport approach HSD, which describes 
successfully both hadronic and electromagnetic observables from p + A and 
A + A collisions from SIS to SPS energies. The production of cc and Drell-Yan 
pairs is treated perturbatively employing experimental cross sections while the 
interactions of cc pairs with hadrons are included by conventional cascade-type 
two-body collisions. Adopting 6 mb for the cc-baryon cross sections the data 
on J/ty suppression in p + A reactions are reproduced in line with calculations 
based on the Glauber model. We study different models for cc dissociation on 
mesons in comparison with the experimental data of the HELIOS-3, NA38 and 
NA50 collaborations. Adopting absorption cross sections with mesons above the 
DD threshold in the order of 1.5 - 3 mb we find that all data on J/ty suppression 
from both proton-nucleus and nucleus-nucleus collisions can be described without 
assuming the formation of a quark-gluon plasma in these collisions. 
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1 Introduction 



The study of hot and dense nuclear matter via relativistic nucleus-nucleus collisions is 
the major aim of high energy heavy-ion physics. Nowadays, the search for a restoration 
of chiral symmetry at high baryon density and temperature or for a phase transition to 
the quark-gluon plasma (QGP) is of specific interest. In this context Matsui and Satz 
H] have proposed that a suppression of the J/^f yield in ultra- relativistic heavy- ion 
collisions is a plausible signature for the formation of the quark-gluon plasma because 
the J/ty should dissolve in the QGP due to color screening |L|. This suggestion has 
stimulated a number of heavy-ion experiments at CERN SPS to measure the J/\& via 
its dimuon decay. Indeed, these experiments have shown a significant reduction of the 
J/ty yield when going from proton-nucleus to nucleus-nucleus collisions 0. Especially 
for Pb + Pb at 160 GeV/A an even more dramatic reduction of J/^> has been reported 
by the NA50 collaboration §, |. 

To interpret the experimental results, various models based on J/ty absorption by 
hadrons have also been proposed. In Ref. ||, Gerschel and Huffier have shown within 
the Glauber model that the observed suppression of J/^f in nuclear collisions is con- 
sistent with the hadronic absorption scenario if one assumes a J/\l/-nucleon absorption 
cross section of about 6-7 mb. Similar but more recent analyses by the NA50 collab- 
oration and Kharzeev || have led to the same conclusion. However, this model 
has failed to explain the "anomalous" suppression reported in central Pb + Pb colli- 
sions, thus leading to the suggestion of a possible formation of a quark-gluon plasma 
in these collisions ||, || 0, @|- On the other hand, Gavin et al. [§, |10| , based also on the 
hadronic absorption model, have found that although J/^f absorption by nucleons is 
sufficient to account for the measured total J/^f cross sections in both proton- nucleus 
and nucleus-nucleus collisions, it cannot explain the transverse energy dependence of 
J/^f suppression in nucleus-nucleus collisions. To account for the nucleus-nucleus data 
they have introduced additionally the absorption on mesons ('comovers') with a cross 



section of about 3 mb. A similar model has also been proposed by Capella et al. [11 



to describe the J/\I/ and ty' suppression in nucleus- nucleus collisions. On the other 



hand, Kharzeev et al. [0 claim the 'comover' absorption model to be inconsistent 
when considering all data on J/^> production simultaneously. 

In all these studies the dynamics of the collisions is based on the Glauber model, 
so a detailed space and time evolution of the colliding system is not included. In 
particular, the transverse expansion of the system and the finite hadron formation time 
is ignored in the Glauber models. Especially for nucleus-nucleus collisions involving 
heavier beams, such as the Pb + Pb collisions at 160 GeV/nucleon, the dynamics 
is more subtle than in proton and S induced reactions. Thus dynamical models are 
needed to complement our information on the reaction dynamics. In this respect Loh 
et al. (TB| have investigated the J/^ dissociation in a color electric flux tube in a 
semiclassical model based on the Friedberg-Lee color dielectric Lagrangian. They find 
that the cc dissociation time is in the order of 1 fm/c. A first transport theoretical 



analysis of J/^f production and absorption has been performed in Ref. |14| where the 



cc production was based on the LUND string formation and fragmentation model 
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Indeed, substantial differences to the Glauber approaches have been found due to a 
finite formation time of the cc pair and its subsequent interactions with baryons and 
mesons. However, due to the low statistics achieved in the numerical calculations a 
definite conclusion about the charmonium suppression could not be obtained in the 
latter study since alternative production and absorption schemes also could lead to 
different absorption rates. 

In this work we continue the studies in Ref. [13] on cc production and suppression 
at SPS energies within the covariant transport approach HSD 1 [[L6[] using different 
production and absorption models. The nonequilibrium transport calculations have 
been shown to describe satisfactorily both the measured hadronic observables (rapidity 
distributions, transverse momentum spectra etc. [16|]), which are sensitive to the final 



stage collision dynamics, and dilepton spectra |T7L p~8L [L9L 20 1, which reflect also the 



initial hot dense stage of the collisions. It thus gives a more realistic description of 
the heavy- ion reaction dynamics than that used in Refs. |^, |], [7], [| M. Within 
this approach we can check if the cc pair might be destroyed by nucleons before the 
mesons are produced as argued in Ref. |7[] (cf. also Refs. || [|). Furthermore, we can 
test if a finite lifetime of the cc pre-resonance system as suggested by Kharzeev et al. 
[[], comes in conflict with the data since according to Ref. |22j the J/\I/ - meson 



cross section might be negligibly small in hadronic matter due to the small size of the 
J/\l/ and its large mass gap from open charms. However, it is expected that the cc pair 
is first produced in a color-octet state together with a gluon ('pre-resonance state') 
and that this more extended configuration has a larger interaction cross section with 
baryons and mesons before the J/\l/ singlet state finally emerges. We will not address 
the question of whether the magnitude of the cc-hadron cross sections used is correct 
or can be justified by nonperturbative QCD and concentrate on the question if specific 
reactions models can be ruled out in comparison to the data available so far. 

The outline of the paper is as follows: In Section 2 we briefly describe the covariant 
transport approach employed and present details of the baryon and meson dynamics for 
S + U at 200 GeV/A and Pb + Pb at 160 GeV/A. Section 3 contains the description 
of charmonium and Drell-Yan production by nucleon-nucleon collisions as well as two 
models for the charmonium reabsorption on hadrons. In Section 4 we present a detailed 
comparison of our calculations with the experimental data available so far while Section 
5 concludes with a summary and discussion of open problems. 



2 The covariant transport approach 



In this work we perform our analysis along the line of the HSD approach which is 
based on a coupled set of covariant transport equations for the phase-space distributions 
fh{ x iV) °f hadron h |16], |23| , i.e. 



{(n, - n u d;u» - Mffiug) d% + (n„a^ + M* h c%u£) d»} f h (x, P ) 



1 Hadron String Dynamics 
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J d2d3d4 . . . [G f G] 12-^34.. £ 4 (n + n 2 - n 3 - n 4 . . .) 

/12A3/14... 

x{fh 3 (x, Pz)fhA ( x , Pi)fh (x, p) f h2 (x,p 2 ) 

- fh ( x ,P)fh 2 (x,p 2 )f h[i ( x >P3) fh 4 ( x >P4)} ■■ • • (!) 

In Eq. ([1]) U^(x,p) and Ufi(x,p) denote the real part of the scalar and vector hadron 
selfenergies, respectively, while 12^34.. .^(Jl + n 2 — n 3 — n 4 . . .) is the 'transition 

rate' for the process 1 + 2^3 + 4 + ... which is taken to be on-shell in the semiclassical 
limit adopted. The hadron quasi-particle properties in (|l|) are defined via the mass-shell 
constraint 

^IF-M* 2 ) , (2) 
with effective masses and momenta given by 

M*(x,p) = M h + U*(x,p) 

U.»(x,p) = i?-Uit{x,p) , (3) 
while the phase-space factors 

f h (x,p) = l±f h (x,p) (4) 

are responsible for fermion Pauli-blocking or Bose enhancement, respectively, depend- 
ing on the type of hadron in the final/initial channel. The dots in Eq. (ffl) stand 
for further contributions to the collision term with more than two hadrons in the fi- 
nal/initial channels. The transport approach ([I]) is fully specified by U^(x,p) and 
U£(x,p) (/i = 0,1,2,3), which determine the mean-field propagation of the hadrons, 
and by the transition rates G^G5 4 (. . .) in the collision term, that describes the scat- 
tering and hadron production/absorption rates. 

The scalar and vector mean fields and Ufc for baryons are taken from Ref. [IB 



and don't have to be specified here again, since variations in the baryon selfenergies 
within the constraints provided by experimental data were found to have no sizeable 
effect on the issue of cc production and absorption. In the present approach we propa- 
gate explicitly - apart from the baryons (cf. ||16|| ) - pions, kaons, 77's, ?/'s, the Y~ vector 
mesons p, u, and K*'s as well as the axial vector meson ai. The production of mesons 
and baryon-antibaryon pairs is treated within the LUND string model employing 
a formation time t> = 0.7 fm/c which controls the baryon and meson rapidity dis- 
tributions dN/dy in comparison to experimental data. As meson- meson channels we 
include the reactions titt — > p, titt —>■ KK, np — > <f), np — > ai as well as the time reversed 



reactions using Breit-Wigner cross sections with parameters from the literature pi 
and exploiting detailed balance. As noted before, this transport approach was found 
to describe reasonably well hadronic as well as dilepton data from SIS to SPS energies 



ig, [19|, |2C 

Before discussing the question of charmonium production we show in Fig. [I] the time 
evolution of the baryon density Pb(x, y, z; t) as a function of z and time t for x — y = 
in a central collision of Pb + Pb at 160 GeV/A in the nucleon-nucleon center-of-mass 
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frame. The two Pb-ions start to overlap at t ~ 1 fm/c, get compressed up to a maximum 
density of about 2.5 fm~ 3 at t ~ 2.5 fm/c and expand in longitudinal (z) direction 
later on indicating a sizeable amount of transparency. We note that the space-time 
evolution in Fig. |l] is controlled by the experimental rapidity distributions dN/ dy for 
protons and negatively charged particles from NA49 |25j ; the streaming of hadrons or 
more precisely their distribution in velocity (3, i.e. dN/d(3 can also be directly extracted 
from the experimental data using dN/d/3 = \dN/dy(j3) (1/(1 + (3) + 1/(1 — (3)). 

The produced mesons in this reaction (in a central cylinder of radius R = 3 f m 
and volume V ~ 15 fm 3 ) appear at about t « 2 fm/c as can be extracted from Fig. [2] 
(lower part) where the densities of pions, p, u and 77 mesons are displayed separately 
as a function of time. The maximum in the meson density (~ 1 fm~ 3 for pions) at t ~ 
3.2 fm/c appears with a delay of tf = 0.7 fm/c with respect to the maximum baryon 
density (cf. Fig. |l|). For comparison we also show in Fig. ^| the meson densities for a 
central S + U collision at 200 GeV/A within the same volume. When comparing to 
the central Pb + Pb collision at 160 GeV/A we observe a lower meson density for the 
S + U case in the central overlap region; especially the p-meson density is lower by 
about a factor of 2. Nevertheless, high baryon and meson densities are encountered in 
these reactions for t > 2 fm/c where a cc pair - that can be produced from 1 - 3.5 fm/c 
- has to pass through. 

This situation is summarized schematically in Fig. |3] for a p + Pb (upper part) and 
S + U collision at 200 GeV (middle part) as well as for a central Pb + Pb collision 
at 160 GeV/A (lower part) for freely streaming baryons (thick lines). In all cases the 
initial string formation space-time points are indicated by the full dots; the mesons 
(indicated by arrows) hadronize after a time delay Tp = 0.7 fm/c as shown by the first 
hyperbola. A cc pair produced in the initial hard nucleon-nucleon collision cannot be 
absorped by mesons in the dark shaded areas in space and time; however, in case of 
Pb + Pb, where cc pairs should be produced within the inner rectangles, a sizeable 
fraction will also be produced in a dense mesonic environment. This fraction of pairs 
produced at finite meson density for S + U is much reduced as can be seen from Fig. ^| 
(middle part). 

The upper hyperbolas in Fig. 3 represent the boundaries for the appearance of 
mesons from the second interaction points (full dots) which appear somewhat later 
in time; they stand for a representative further nucleon - nucleon collision during the 
reaction. 



3 Charmonium and Drell-Yan production 

Since the probability of producing a cc or Drell-Yan pair is very small, a perturbative 
approach is used for technical reasons. Whenever two nucleons collide a cc pair is 
produced with a probability factor W, which is given by the ratio of the J/^> (or ty') 
to NN cross section at a center-of-mass energy y/s of the baryon-baryon collision, 

W = ff **^/*+*(Vf) (5 ) 

&BB->BB+x{\/ S) 
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The parametrization used for the J/\I/ cross section is 



Vbb^j/*+x{V~s) = 160 ( X " ^) ^ H ( 6 ) 

with a = 7.47 GeV which is shown in Fig. f| (solid line) in comparison to the experi- 
mental data displayed in Ref . [EB| . For the systems to be studied, in the energy regime 



y/s < 30 GeV, the estimated uncertainty of our parametrization is about 20 %. We 
note that the parametrization used conventionally |2T| is 
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^pp^j/'f+xiVs) = d \ l- -j=j (7) 

with c = 3.097 GeV, d = 2 ■ 37 /B^ nb while = 0.0597 is the branching ratio of the 
J/\l/ to dimuons. The expression (|7|) (dashed line in Fig. £|) is in a good agreement with 
our formula (|6]) for 10 GeV < \/s < 30 GeV which is the energy regime of interest. 
The rapidity distribution of the cc-pair is approximated by a Gaussian in the nucleon- 
nucleon center-of-mass of width a ~ 0.6 as in Ref. |27fl while the transverse momentum 



distribution is fitted to experimental data (see below). For production we employ 
the same model, however, scale the experimental cross section by a factor of 0.122 
relative to the J/\l/ production cross section. 

In extension to Ref. |14| the Drell-Yan process is taken into account explicitly. The 



generation of Drell-Yan events was performed with the PYTHIA event generator 



version 5.7 using GRV LO |29) or MRS A [[KJ structure functions from the PDFLIB 
package |JT[] with = 1.0 GeV. In the kinematical domain 1.5 < M < 5.5 GeV, 
3 < yiab < 4 and —0.5 < coshes < 0.5, where 6cs is the polar angle in the Collins- 
Soper reference frame, this yields a dimuon cross section of 270 pb in pp collisions at 
200 GeV and 261 pb in pn collisions, respectively, as in Refs. 0, j27f . 

According to our dynamical prescription the Drell-Yan pairs can be created in each 
hard pp,pn, np or nn collision (y'i > 10 GeV). Since PYTHIA calculates the Drell-Yan 
process in leading order only (using GRV LO structure functions) we have multiplied 
the Drell-Yan yield for iVW collisions by a K-factor of 2.0 (cf. Refs. 0, p7|). In case 
of MRS A structure functions, which include NLO corrections, a K-factor of about 1.6 
had to be introduced (cf. |26[]). The energy distribution of the hard iViV collisions 
dN/dy/s for nucleus- nucleus collisions in our transport approach shows a pronounced 



peak around y/s^ = y2m N (T kin + 2m^) where T kin is the kinetic energy per nucleon 
in the laboratory frame. Thus the main contribution to the dimuon yield - summed 
over all NN events - comes from NN collisions with y/s ~ y/so- However, there are 
also Drell-Yan pairs from iViV collisions with \fs larger or smaller than y/so. 

We have compared our results within the production scheme given above with that 
used in Refs. 0, |3|, |27] , where the Drell-Yan yield from p + A and A + A collisions is 
calculated as the isotopical combination of the yield from pp and pn at fixed ^/so scaled 
by A P x A T . We found that the variation from the scheme used in Refs. 0, [3], 



is less then 10%. We, furthermore, note that the difference in the dimuon spectrum 
using different structure functions (GRV LO or MRS A) is less then 5%. In the present 
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analysis we have discarded dimuon production from open charm channels because a 



recent analysis by Braun-Munzinger et al. [32] on the basis of the same PYTHIA event 



generator has shown that the open charm contributions at low and high invariant 
masses are of minor importance. 

Since the production scheme is the same for cc-pairsQ their total cross section (with- 
out reabsorption) also scales with Ap x At', the ratio of the J/ty to the Drell-Yan cross 
section thus provides a direct measure for the J/\l/ suppression. 

In order to obtain some information about the primary distribution of the produced 
'pre-resonance' states in coordinate space we show in Fig. |5| the cc distribution in the 
(x, z) -plane integrated over y for a central collision of Pb + Pb at 160 GeV/A. The 
collision of the two nuclei proceeds along the ^-direction and the actual z-axis has been 
stretched by the Lorentz factor 7 cm ^ 9.3 to compensate for the Lorentz contraction in 
beam direction. It is clearly seen from Fig. § that the production of the charmonium 
state by hard nucleon-nucleon collisions is enhanced in the center and drops rapidly in 
the surface region of the overlapping nuclei. 

We follow the motion of the cc pair in hadronic matter throughout the collision dy- 
namics by propagating it as a free particle. In our simulations the cc pair, furthermore, 
may be destroyed in collisions with hadrons using the minimum distance concept as 
described in Sec. 2.3 of Ref. For the actual cross sections employed we study two 
models (denoted by I and II) which both assume that the cc pair initially is produced 
in a color-octet state and immediately picks up a soft gluon to form a color neutral 
cc — g Fock state (color dipole). This extended configuration in space is assumed 
to have a 6 mb dissociation cross section in collisions with baryons (cc + B —>■ A c + D) 
as in Refs. |4], [5], § during the lifetime r of the cc — g state which is a parameter. In 
the model I we assume r = 10 fm/c which is large compared to the nucleus-nucleus 
reaction time such that the final resonance states J/\l/ and ty' are formed in the vac- 
uum without further interactions with hadrons. In the model II we adopt r = 0.3 
fm/c as suggested by Kharzeev |J which implies also to specify the dissociation cross 
sections of the formed resonances J/ty and Nl/' on baryons. For simplicity we use 3 mb 
following Ref. [[Hp. The cross section for cc — g,J/^, or dissociation on mesons 
(cc + m — > DD) is treated as a free parameter ranging from to 3 mb. 

In order to test the transverse momentum dependence of the J/^f production we 
have performed calculations for p + U and S + U at 200 GeV/A (using model I 
with o h abs° ns = 6 mb and <J™ s sons = 1 mb) which are scaled in magnitude to the 
respective data for the same systems from Ref. 0. Our calculated results are shown 
in Fig. |G] for both systems in terms of the histograms while the solid lines represent fits 
to the experimental transverse momentum distributions. Since this pt dependence is 
described fairly well we expect the cc event distributions to be quite realistic. 



4 Analysis of experimental data 



2 The formation time tf for the cc-pair is assumed to be zero in our present approach in contrast 
to Ref. @. 
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4.1 HELIOS-3 data 



Since the J/\I/ and ^' are measured in nuclear collisions through their decay into 
dimuons, we calculate explicitly the dimuon invariant mass spectra from these colli- 
sions. This includes not only the decay of the J/^,^' but also the decay of other 
mesons (77, p, u>, and 0) as well as the Dalitz decay of the 77, u>, rf etc. Details on cal- 
culating the dimuon spectra from heavy-ion collisions up to invariant masses of about 
1.3 GeV can be found in Refs. |17|, [18], |(J. Since the Drell-Yan contribution is im- 
portant for dileptons with invariant masses above 1.5 GeV pTJ we have included these 



channels by computing for each 'hard' nucleon-nucleon collision > 10 GeV) their 



contribution via PYTHIA 5.7 |28| as described above. 

We have carried out calculations for p + W and central S + W collisions at 
200 GeV/A as in Ref. |TJ| in order to check on an absolute scale if the production and 
absorption schemes employed in the HSD approach are included properly. In Fig. [7| 
we show the dilepton invariant mass spectra for p + W and S + W at 200 GeV/A 
normalized to the number of charged particles in the pseudorapidity bin 3.7 < rj < 5.2 



and compare them with the experimental data from Ref. [34| including their accep 



tance and mass resolution. Again we have employed model I with cr^l r s yons = 6 mb and 
<T X SC " 1S = 1 m b ( see below). It is seen from Fig. |7| that in the p + W and S + W cases 
the theoretical results agree well with the data on an absolute scale which implies that 
apart from the low mass dimuon spectrum - which has been analyzed in Ref. - also 
the J/\l/ and region is described reasonably well. At low dimuon masses the explicit 
contributions from the mesons rj, u>, <fr are displayed in terms of the thin lines while for 
invariant masses above 1.5 GeV the Drell-Yan (DY), J/^> and fy' contributions are 
shown explicitly. The full solid curve represents the sum of all contributions (without 
open charm channels). It is intersected from M « 1.2 GeV to 1.5 GeV because the 
continuation of the Drell-Yan contribution to lower energies is not clear as well as the 
tail of the p-meson contribution at higher M. Futhermore, open charm channels (from 
D and D mesons) have not been included yet (cf. Ref. [0). We note that for S + W 
in the invariant mass range from 1.6 GeV < M < 2.5 GeV the dimuon yield is al- 
most twice that for p + W within the present normalization and that our calculations 
reproduce also the intermediate mass range on the basis of PYTHIA 5.7 reasonably 
well. Thus our calculations can be regarded as fully microscopic studies on dilepton 
production for invariant masses 0.2 GeV < M < 5 GeV including both the relevant 
'soft' and 'hard' processes. 

4.2 Proton-nucleus reactions 

Since the absorption of cc-pairs on secondary mesons in proton-nucleus collisions is 



practically negligible these reactions allow to fix 'experimentally' the cc-baryon 
dissociation cross section. Furthermore, the Glauber approach for cc absorption in 
this case should be approximately valid at energies of about 200 - 450 GeV such that 
the transport calculations can be tested additionally in comparison e.g. to the model 
calculations of Refs. §, |, % |TU|, [TTJ. 
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In Fig. [8] we show the results of our calculations for the J/^f 'survival probability' 
(open squares) using 6 mb for the absorption cross section of the cc-pairs on nucleons 
in comparison to the data ||. The experimental 'survival probabilities' S^f in this 
figure as well as in the following comparisons are defined by the ratio of experimental 
J/^J to Drell-Yan cross sections as 

/ R J/% 

K&AB 1 2.9-4.5 GeV 

where A and B denote the target and projectile mass while a J J^ and <j1b stand for 
the J/ty and Drell-Yan cross sections from AB collisions, respectively, and B w is the 
branching ratio of J/^> to dimuons. The theoretical ratio is defined as 

°theor ~ a j > \ J ) 

where Nj/y is the multiplicity of initially produced J/\I/'s while Mj/q, is the multiplicity 
of J/^'s that survive the hadronic final state interactions. We note that due to the 
large statistical error bars of the experimental data absorption cross sections <J b ^ s yons of 
6 ± 1 mb are compatible also. These values are slightly smaller than those of Kharzeev 
et al. |T^] in the Glauber model claiming 7.3 ± 0.6 mb, but in the same range as 



those used in the Glauber models of Ref. ||. We do not expect to get exactly the 
same values as in the Glauber model because the transverse expansion of the scattered 
nucleons as well as of the cc-pairs is neglected there. Due to the 'optimal' reproduction 
of the available data with a cross section of about 6 mb we will use this value also for 
nucleus-nucleus collisions in the following. 

4.3 Transverse energy distributions 

In order to perform a detailed comparison to the data of the NA38 and NA50 Collab- 
orations for S + U at 200 GeV/A and Pb + Pb at 160 GeV/ A we first have to fix the 
experimental event classes as a function of the neutral transverse energy Et in order to 
allow for an event by event analysis. In this respect we compute the differential cross 
section 

da r°° , „ dN 



2tt / b db-—(b) (10) 
dE T Jo dE T ' v ; 

as a function of the impact parameter b, where ^^(&) is the differential Et distribution 
for fixed b. Since the detector response is not known to the authors the Et distribution 
is rescaled to reproduce the tail in the experimental Et distributions which results in 
the dashed histograms in Fig. |] that increase for low Et- 

In the actual experiments, however, only events are recorded with a yU + /i~ pair of 
invariant mass M > 1.5 GeV and rapidity 3 < y iab < 4 for NA38 or M > 2.9 GeV 
and rapidity 2.93 < yi a b < 3.93 for NA50. A respective selection in the transport 
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calculation is obtained by 



- 2vriV y bdb YlWr(b), (11) 
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where W^{b) are the weights for produced p/i pairs within the experimental cuts. iVo is 
a normalization factor to adjust to the experimental number of events. The calculated 
distributions (|ll|) are shown in Fig. |9| in terms of the solid histograms which reasonable 
reproduce the experimental distributions (grey histogram for S + U, full dots for Pb 
+ Pb). The average values for Et in the 5 bins are indicated in Fig. |9] by the open 
squares and coincide with the corresponding experimental values from Ref. [[|. 



4.4 NA38 and NA50 data 

We are now in the position to perform a comparison to the experimental survival 
probabilities for J/\l/ and production in the respective transverse energy bins. We 
first compute the results for S + U at 200 GeV/A and Pb + Pb at 160 GeV/A within 
the model I varying the dissociation cross section on mesons of the cc — g object from 
to 1.5 mb while keeping the absorption cross section on baryons fixed at 6 mb. The 
calculated Sj[f or are displayed in Fig. |10| in comparison to the data for both systems; 
the dashed lines are obtained for cr™ s sons = mb while the solid lines correspond to 
<J™ s sons = 1.5 mb. Whereas the data for S + U appear to be approximately compatible 
with our calculations without any dissociation by mesonsQ the Pb + Pb system shows 
an additional suppression. This finding is in agreement with the results of Glauber 
models |j[ [T2j . On the other hand, the Pb + Pb data are well reproduced with a cross 
section of 1.5 mb (in model I) for the J/\l/ absorption on mesons which, however, then 
slightly overestimates the suppression for the S + U data for the 3 middle Ej- bins. 
Our calculations thus do not indicate a strong argument in favor of a QGP phase in 
the Pb + Pb reaction to interpret the J/\l/ suppression within the model I. 

The latter conclusion is different from the Glauber calculations of Refs. [0, [12|] and 
should be due to the simplified assumptions about the actual meson abundancy with 
which the cc pair can be dissociated. Note again that mesons appear in our dynamical 
approach only after Tp = 0.7 fm/c after the first 'hard' collision. As has been pointed 
out by Gavin et al. M especially the J/\l/ in the comoving frame should be dissociated 
by p and u mesons and less by pions due to the large gap in energy for DD dissocia- 
tion. We have investigated this suggestion more quantitatively within our microscopic 
approach and find the following hadronic decomposition for J/ty absorption on mesons: 
S + U (central): pions (35%), p's (42%), c/s (15%), 77's (8%); Pb + Pb (central): 
pions (37%), p's (42%), cj's (13%), 77's (8%). Thus the hadronic decomposition prac- 
tically does not change when going from S + U to the heavier Pb + Pb system. As a 
side remark we note that the dilepton yield from central S + Au and Pb + Au collisions 
- normalized to the charged particle multiplicity in the rapidity bin 2.1<y<3.1- 
experimentally appears to be the same which is in line with our findings. 



3 The S+U data are best described with cj'°.^ vons w 6.5 mb. 
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A remarkable difference, however, is found when comparing the number of absorbed 
J/\l/'s by mesons in central collisions for S + U and Pb + Pb as a function of time. 
In order to compare the two absorption rates dN^ s (t)/dt on a relative scale we have 
multiplied the absorption rate for the system S + U by the ratio of projectile nucleons 
(208/32) in Fig. [TI[ Here the heavier system Pb + Pb shows higher absorption rates 
from the beginning (at t = 2 fm/c), which correlates with the central meson densities 
shown in Fig. and also the J/ty absorption on mesons lasts longer in accordance with 
the calculations in Ref. fL4| . The latter effect can also be extracted from the schematic 
picture in Fig. |3] where the number of produced cc-pairs at finite meson density is 
expected to be much larger for Pb + Pb than for S + U. 

We, furthermore, compute the ratio S J ^ for S + U at 200 GeV/A and Pb + Pb 
at 160 GeV/A within the model II for a cc — g lifetime r = 0.3 fm/c varying the 
dissociation cross section of the cc— pair with mesons from 0-3 mb while keeping 
the absorption cross section on baryons fixed at <7^ yons = 6 mb for the 'pre- resonance' 
state and at 3 mb for the formed J/^f resonance. The calculated ratios are displayed in 
Fig. [12| in comparison to the data; again the dashed lines correspond to the calculations 
without any charmonium absorption on mesons whereas the solid lines represent our 
calculations for a meson absorption cross section of 3 mb. In the absorption model II the 
data for S + U are no longer compatible with our calculations without any dissociation 
by mesons. The S + U data here need an absorption by mesons in the range of 3 mb 
as in the phenomenological model of Gavin et al. || |10[ . With (T^ s sons ~ 3 mb for the 
absorption on mesons, however, the Pb + Pb data appear to be compatible, too. For 
future experimental studies we display in Fig. [12] the calculated ratio Sf[^ or (within 
model II) for more peripheral reactions (open diamonds) that correspond to an E^-bin 
from 5^10 GeV in case of S + U at 200 GeV/A and to an £ T -bin from 10 ^ 20 GeV 
in case of Pb + Pb at 160 GeV/A. 

We additionally investigate if the ^' suppression measured by NA38 and NA50 for 
S + U and Pb + Pb can be described simultaneously within our approach without 
including any additional assumptions. Since the suppression in proton-nucleus re- 
actions is practically the same [2fJ we employ also a dissociation cross section of 6 mb 
on baryons while the absorption cross section on mesons is treated again as a free pa- 
rameter. Our numerical results for both systems on the ty' survival probability within 
the absorption model I are displayed in Fig. [13] in comparison with the data from Q . 
Due to the higher suppression in these reactions the absorption on mesons requires 
cross sections from 2.5 - 3.5 mb. We note that we cannot describe the experimental 
ratios Sf xp for central S + U collisions (high Et) within the present treatment which 
might be due to the neglect of — > J/^> + X, or tt^' — > J/ty + X reaction channels. 
This also holds for calculations within the model II where cr^l e s sons ~ 4 — 5 mb has to be 
assumed. Furthermore, the actual statistics reached for production and absorption 
events is too low to allow for final conclusions here. 



11 



5 Summary 



In this paper we have carried out a microscopic transport study of J/^f, and Drell- 
Yan production in proton-nucleus and nucleus-nucleus collisions. Our calculations show 
that the absorption of 'pre-resonance' cc — g states by both nucleons and produced 
mesons can explain reasonably not only the inclusive J/ty cross sections but also the 
transverse energy (Et) dependence of J/ty suppression measured in nucleus-nucleus 
collisions. In particular, the absorption of J/^'s by produced mesons is found to be 
important especially for Pb + Pb reactions, where the J/\l/-hadron reactions extend 
to longer times as compared to the S + W or S + U reactions. This is in contrast 
with results based on a simple Glauber model, which neglects both the transverse 
expansion of the hadronic system and the finite meson formation times, where the 
cc — N absorption is roughly sufficient even for S- induced collisions. As a consequence 
we do not find a necessary argument to require the formation of a quark-gluon-plasma 
in Pb + Pb collisions at 160 GeV/A. This could only be done through experimental 
or theoretical evidence that the charmonium-meson cross sections employed (« 1.5 - 
3 mb) in our analysis are too large. Since the cc dissociation on mesons is expected to 
be dominated by flavor exchange reactions the present cross sections in our opinion, 
however, should be reasonable. 

We close our study by noting that in the phase of high baryon and meson density 
both the meson densities (cf. Fig. ||) as well as the associated energy densities |14| are 
very high such that a purely hadronic reaction scheme might be questionable. Further 
experimental data with good statistics also for p + A collisions with light nuclei are 
expected to provide more accurate constraints. 
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Figure 1: The baryon density ps{x = 0,y = 0,z;t) for a Pb + Pb collision at 
160 GeV/A and impact parameter 6 = 1 fm. 
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Figure 2: The density of ir, p, uj and r] mesons in a central cylinder of radius R = 3 fm 
and volume V ~ 15 fm 3 for a S + U collision at 200 GeV/A and a Pb + Pb collision 
at 160 GeV/A for b = 1 fm. 
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Figure 3: Schematic representation of p + Pb collision at 200 GeV (upper part), S + U 
collision at 200 GeV/A (middle part) and a Pb + Pb collision at 160 GeV/A (lower 
part) in space-time. The full dots represent early hard collision events (for Drell-Yan 
and cc pairs) while mesons (7r,rj,p,uj, etc. - arrows) only appear after a respective 
formation time Tp ~ 0.7 fm/c. The overlap area (inner rectangle) specifies the space- 
time region of hard production events. 
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Figure 4: The parameterization (|6]) for the elementary J/ty cross section in pp collisions 
(solid line) in comparison to the experimental data from Ref. |26|]. The dashed line 
represents the parametrization (|7|). 
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Figure 5: The production probability for cc-pairs in a central Pb + Pb collision at 
160 GeV/A in the (x, z)-plane integrated over y. The z-axis is scaled by the Lorentz 
factor 7 cm = 9.3 to compensate for the Lorentz contraction in beam direction. 
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Figure 6: The transverse momentum distribution of J/^f mesons in p + U (upper 
part) and S + U collisions (lower part) at 200 GeV/A from the HSD calculations 
(histograms). The solid lines are fits to the experimental data from Ref. @. 
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Figure 7: Dimuon invariant mass spectra from p + W and S + W collisions at 
200 GeV/nucleon in comparison to the data of the HELIOS-3 collaboration P^| . In the 
low mass region the individual contributions from the rj Dalitz decay, u and $ decay 
are shown by the thin lines. The Drell-Yan contribution (denoted by DY) is calculated 
only for invariant masses M > 1.5 GeV. The explicit contributions from J/ty and \P' 
decays are indicated by thin lines for M > 2.7 GeV. The thick solid lines represent the 
sum of all dimuon channels (except for open charm contributions). 
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Figure 8: The J/\& 'survival probability' S J ^ for p + A reactions at 200 GeV assuming 
a 6 mb cross section for the cc dissociation on baryons (open squares) in comparison 
to the experimental data from within the model I. 
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Figure 9: The distributions in the neutral transverse energy for S + U at 200 GeV/A 
and Pb + Pb at 160 GeV/A; dashed histograms: result of the transport calculations 
without any constraints on the experimental acceptance (10); experimental distribu- 
tions: grey histogram for S + U, full dots for Pb + Pb; solid histograms: HSD cal- 
culation for the Et distribution according to Eq. (|TT|). The open squares denote the 
average Et in the experimental Et bins. 
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Figure 10: The ratio S J/ * for S + U at 200 GeV/ A (upper part) and Pb + Pb at 
160 GeV/A (lower part) as a function of the transverse energy in comparison to the 
experimental data from within the model I assuming a long lifetime for the cc — g 
system. The absorption cross section on mesons is varied from (dashed lines) to 
1.5 mb (solid lines). 



24 




Figure 11: The J/\l/ absorption rate dN^ s (t)/dt on mesons for central collisions of 
S + U at 200 GeV/A and Pb + Pb at 160 GeV/A. The absorption rate for S + U has 
been multiplied by a factor 208/32 to compensate for the different number of projectile 
nucleons. 



25 



1.0 



0.8 



0.6 



model II 



8: 

i-s 



0.4 



0.2 



0.0 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



1 1 1 1 

- 


S+U 200 GeV/ A 






# 




- 


• — t — . 


a abs ba ^" s (pre-res.) = 6mb 




a ab s bary0nS CJ/"*) = 3 mb 




a . mesons = - 3 mb 

abs 

i , i 


i , i 



20 



40 60 

E T (GeV) 



80 



Pb+Pb, 160 GeV/A 



o 



tf abs baryons (pre-res.) = 6 mb 

a ab s bary ° nS CJ'*) = 3 mb 
a „ mesons = - 3 mb 



20 40 60 80 100 

E T (GeV) 



120 



140 160 



Figure 12: The ratio S J/qi for S + U at 200 GeV/A (upper part) and Pb + Pb at 
160 GeV/A (lower part) as a function of the transverse energy in comparison to the 
experimental data from || within the model II (see text). The absorption cross section 
on mesons is varied from (dashed lines) to 3 mb (solid lines); the dissociation cross 
section on baryons for the pre-resonance cc — g system was taken as 6 mb while for the 
J/^f singlett cross section with baryons 3 mb were adopted. The open diamonds repre- 
sent the calculated 'survival probabilities' S J '^ (within model II) for more peripheral 
reactions (see text). 
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Figure 13: The ratio S*' for S + U at 200 GeV/ A (upper part) and Pb + Pb at 
160 GeV/A (lower part) as a function of the transverse energy in comparison to the 
experimental data from within the model I assuming a long lifetime for the cc — g 
system. The absorption cross section on mesons is varied from 2.5 to 3.5 mb. 
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